Tailoring of vibrational state populations with light-induced potentials in molecules 
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We propose a method for achieving highly efficient transfer between the vibrational states in a 
diatomic molecule. The process is mediated by strong laser pulses and can be understood in terms 
of light-induced potentials, fn addition to describing a specific molecular system, our results show 
how, in general, one can manipulate the populations of the different quantum states in double well 
systems. 
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Quantum control of molecular processes has recently 
opened new possibilities for controlling chemical pro- 
cesses, and also for understanding the multistate quan- 
tum dynamics of molecular systems. With strong and 
short laser pulses one can create quantum superpositions 
of vibrational and continuum states, i.e., wave packets, 
which often propagate like classical objects under the in- 
fluence of molecular electronic potentials 0]. Typically 
in spectroscopy one induces with laser light population 
transfer between individual vibrational states, and de- 
scribes the process using Franck-Condon factors, pertur- 
bation theory and a very truncated Hilbert space for the 
molecular wave functions. If the laser light comes in as 
a very short pulse we expect to couple many vibrational 
states because of the broad spectrum. However, in this 
Letter we show that it is possible to start from a single 
vibrational state, and end up selectively on another sin- 
gle vibrational state, with very high efficiency, even when 
one uses strong and fast laser pulses. 

We have previously discussed how one can transfer 
the vibrational ground state population of one electronic 
state to the vibrational ground state of another elec- 
tronic state [pi. This process was called adiabatic pas- 
sage in light-induced potentials (APLIP). These poten- 
tials, which depend on time, due to the time dependence 
of the laser pulse envelopes, provide a useful description 
of the process (see e.g. Ref. || and references therein). 
Here we use the same kind of light-induced potentials to 
describe and understand transfer processes between ex- 
cited vibrational states. However, we will see that the 
process is more complicated than the simple adiabatic 
following assumed in APLIP. 

To demonstrate the new process we have chosen the 
same Na dimer potentials, shown in Fig. |l|(a), as in our 
previous study 0. The three electronic states are cou- 
pled in a ladder formation by two laser pulses. Instead of 
describing the system in terms of the vibrational states 
and Franck-Condon factors, we consider only the elec- 
tronic state wave functions, ^i(R,t), i — X, A, H. If wc 
apply the rotating wave approximation, we can 'shift' the 
X 1 E+ and 2 1 n s state potentials in energy by the corre- 
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FIG. 1. (a) The three Na2 potential energy surfaces used 
in our calculations: X X E+, A X E+ and 2 1 IT 9 . (b) The shifted 
potentials. Here we have set A = Ai = — A2, [see Eq. (Q)] 
and show the position of the intermediate A X E+ state (dotted 
line) in respect to the other two electronic states (solid lines) 
for HA = ±2000 cm" 1 . 



sponding laser photons. Then we obtain the situation 
shown in Fig. [j](b) , after a suitable redefinition of the en- 
ergy zero point. The resonances between the electronic 
state potentials now become curve crossings. The evolu- 
tion of the three wave functions *i>i(R,t) is given by the 
time-dependent Schrodinger equation with the Hamilto- 
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where R is the internuclear separation, m is the reduced 
mass of the molecule, and the electronic potentials and 
couplings are given by 
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U{R,t) = 



U X (R) hfli(t) 

mi(t) u A (R) + hA 1 
o nn 2 (t) 
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U n {R) + fi(Ai + A 2 ) 
(2) 



Here Ux(R), Ua{R), and Un(R) are the three potentials, 
Ai and A 2 are the detunings of the two pulses from the 
lowest points of the potentials [dashed lines in Fig. 0(a)] , 
and fii(t) = fixAEi(t)/h, fl 2 (t) = ^AivE 2 (t)/fi are the 
two Rabi frequencies. We have assumed for simplicity 
that the two dipole moments are independent of R and we 
have used Gaussian pulse shapes, fij(i) = S!exp{ — [(t — 
U)/T} 2 }, i = l,2. 

The light-induced potentials are obtained by diagonal- 
ising the potential term (0). For these potentials the 
curve crossings become avoided crossings. It is easy to 
see from Fig. 0(b) that in the absence of the pulses (or 
when they both are weak) one of the light-induced states 
corresponds, at low energies, to the double well structure 
formed by the X l H+ and 2 1 II ff electronic state poten- 
tials. In Fig. 0(b) this corresponds to the lower part 
of the two solid curves and we will call this the active 
eigenstate. In Ref. [0 we showed that, if the pulses are 
applied in a counterintuitive order (t\ > t 2 ), the double 
well structure will disappear as the bottom of the (ini- 
tially empty) well on the right moves up and vanishes. 
(In Fig. 0(b) this is seen to be because of a repulsion 
between the 2 1 II g and A 1 ^ (A < 0) state). After this 
the left well broadens and moves to the right. Finally 
the double well structure is re-established as the pulses 
reduce in intensity. If we now consider that the ground 
vibrational state of the left well is initially populated, 
that population would follow the light-induced potential 
and be transformed smoothly into the ground state wave 
function of the right well. This is APLIP, and the smooth 
change in the total wave function is demonstrated in the 
contour plot in Fig. 0(a). 

If we choose Ai = — A 2 = A, the bottoms of the two 
wells are on the same level initially and finally [as in 
Fig. 0(b)]. The sign of A determines the evolution of the 
light-induced potential. For negative A we obtain the 
APLIP situation, but for positive A the right well drops 
down at first, instead of disappearing. In this case the 
APLIP situation is not obtained. 

One is tempted to associate the APLIP process with 
STIRAP and consider the active light-induced po- 
tential as a dark state |J , which, because of the counter- 
intuitive pulse order, remains uncoupled from the other 
states. It is true that during the APLIP process the in- 
termediate state (A 1 E+) population remains very small. 
However, the process is not the same as STIRAP, because 
the strength of the atom-light coupling, i.e. Rabi fre- 
quency, becomes much larger than the vibrational spac- 
ing; it is not possible to isolate only a few energy levels. 
Furthermore, in a STIRAP process there would not be 
a smooth transport of the wave packet in Fig. 0(a) be- 
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FIG. 2. (a) The APLIP process for hA = 
HQ = 2000 cm -1 , h = 25.5 ps, t 2 = 20.5 ps 
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= -3000 cm" 1 , 
and T = 5.5 

ps. (b) The active light-induced state (independent of A at 
low energies) when the pulses are absent (SI = 0). We also 
show the seven lowest vibrational eigenstates \4>i{R)\ 2 of the 
potential. The left well corresponds to X 1 E+ and the right 
one to 2 1 n o . 



cause only two simple vibrational eigenfunctions would 
be involved. 

STIRAP is also considered to be an adiabatic process, 
which is not strictly the case for Fig. 0. We can see this 
by considering the initial (and final) vibrational states of 
the active light-induced potential which matches the two 
original electronic potentials as shown in Fig. 0(b). The 
lowest states correspond to the individual vibrational 
states of the original potentials. In STIRAP we should 
start with the lowest state on left, and we would expect it 
to evolve adiabatically into the lowest state on the right, 
though not in the smooth way seen for the APLIP process 
in Fig. 0(a). However, if the evolution were truly adia- 
batic the population would remain in the left well; by 
definition of adiabatic following it cannot change state. 
This is because, for the example shown in Fig. 0(b), as the 
right well rises up, and the right-hand-well wave function 
with it, there is a diabatic crossing of the energy states 
which allows the population to remain in the left well as 
the right well vanishes. As long as the barrier between 
the two wells is thick enough, the crossing of vibrational 
eigenenergies remains diabatic (uncoupled) and the wave 
packet (left-hand well vibrational ground state) can be 
channelled from left to right (APLIP). There always has 
to be at least one such diabatic crossing, but if there is 
only one, it could come at the start, or the end, of the 
time evolution depending on whether the initial or final 
vibrational state is lower in energy. Thus the interesting 
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processes mediated by the light-induced potentials are 
different from STIRAP, and they can not be obtained by 
perfect adiabatic following. 

Figure ||(b) also displays the initial and final situ- 
ation if we now switch to a positive detuning. Wc 
choose TiA = 2010 cm -1 for the time evolution shown in 
Fig. ||(a) where fiVL = 733 cm -1 and t\ — t<i = T = 5.5 ps. 
We denote the vibrational quantum number of the states 
in the electronic potential with v, and indicate with the 
terms "left well" and "right well" , whether these states 
correspond to the initial state (X 1 S+), or the final state 
(2 1 n g ), respectively. We transfer the v = state on the 
left well into the v = 1 state on the right well. The pro- 
cess is very efficient: the occupation probability of the 
right well (2 1 IT g state) in the end is 94 %. Clearly the 
crucial moment is around t — 20 ps where the original 
single-peak distribution extends to the right and forms 
temporarily a four-peaked distribution over the combined 
well system. Also, around t = 25 ps the four-peak struc- 
ture compresses into the two-peak structure located on 
the right well. 
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FIG. 3. (a) The vibrational state tailoring with 

HA = 2010 cm" 1 , HQ. = 733 cm" 1 , ti = 25.5 ps, t 2 = 20.5 
ps, and T = 5.5 ps. Here we show the time evolution of 
the total probability distribution *(i?,t)| 2 = JT *;(7?,t)| 2 . 
(b) The time-dependent eigenenergies of the five lowest vi- 
brational states of the active light-induced potential for the 
situation shown in (a). [Note that for t < 15 ps and t > 35 
ps the eigenstates match those seen in Fig. ^|(b).] 

In order to understand the process, we have calculated 
the time dependence of the eigenenergies of the active 
light-induced potential [seen for initial and final time in 
Fig. |^(b)] . The figure shows clearly that if the evolution is 
fully adiabatic, no change of v can occur during the time 
evolution because there are no crossings of eigenenergies. 
This also means that no change of well is possible either. 



However, we do see several avoided crossings between the 
eigenenergies and the point is that when the pulses are 
weak, these crossings are passed diabatically, and when 
the pulses are strong they are passed adiabatically. 

In our example, Fig. ||(b), the second lowest state cor- 
responds initially, and finally, to the lowest vibrational 
state of the left well (i.e., the state populated initially 
in our calculation) as may be confirmed by inspecting 
Fig. ^(b). During the time evolution the wave function 
passes the two first avoided crossings diabatically, but for 
t > 20 ps it follows nearly adiabatically the fourth eigen- 
state. Just before t = 30 ps it moves diabatically onto the 
third eigenstate, which asymptotically corresponds with 
the v = 1 state on the right well [see again Fig. §(b)]. 
In Fig. U we show the change at the crucial moment 
around t = 21 ps. The transition from a narrow single 
well wave function into a wide combined two-well wave 
function takes place very quickly. A fascinating aspect 
is that the evolution of the total wave function still fol- 
lows the fourth vibrational eigenstate of the active light- 
induced potential. A similar adiabatic following happens 
also near t = 25 ps when the third and the fourth eigen- 
state again approach each other. 

t = 20.5ps t = 21.0ps t = 21.5ps 
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FIG. 4. The behaviour of the eigenstates of the active 
light-induced potential in the situation of Fig. |^. We see how 
the fourth eigenstate deforms within 1 ps from the initial sin- 
gle-peak structure into the four-peak one. An inverse process 
happens to the third eigenstate. In the diabatic situation the 
system would follow the single-well eigenstates, thus jumping 
from the fourth state to the third. However, here the system 
instead follows the fourth state nearly adiabatically. 

The example given here is not unique. By allowing 
Ai ^ — A2 we have a very good handle in choosing which 
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single well states are paired into a combined well state. 
In Fig. |5| we show two other examples. They demonstrate 
that the process does not require the initial state on the 
left to be v = 0, and the process can be used also to 
change v but keep the initial and final well the same |j| . 
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FIG. 5. The state tailoring process with v — 1 in the left 
well as the initial state, (a) v = 2 in the right well as the 
final state. Here 71 Ai = 2010 cm" 1 and hA 2 = -1960 cm" 1 , 
(b) v — 2 in the left well as the final state. Here hAi = 2010 
cm" 1 and T1A2 = —2160 cm" 1 . 

Even in the region where both pulses are strong the 
adiabatic following of a single vibrational eigenstate is 
not complete. Small oscillating contributions appear in 
the final state. The contour plots tend to emphasise 
these oscillations. For example, in the particular case of 
Fig. ||(a) , only about 2 % probability of diabatic transfer 
creates these oscillations. 

The examples shown here for tailoring the vibrational 
state populations have two advantages over the APLIP 
process. Firstly, having A > (a region in which APLIP 
does not work) moves the central frequencies of the two 
laser pulses away from each other. Secondly, the required 
intensities are clearly smaller as we mainly need to shift 
the potentials, but not to deform them strongly. Of 
course, the discussion in Ref. || regarding the role of 
electronic states outside our truncated three-state sys- 
tem, and the rotational states still apply. We have used 
Gaussian pulse shapes, but with other pulse shapes it 
might be possible to achieve even better control of the 
adiabaticity at different avoided crossings between the 
eigenstates. 

In this Letter we have shown how, even with strong 
and short laser pulses one can still perform very selec- 



tive and yet efficient tailoring of vibrational state pop- 
ulations in molecules. If we were to describe the pro- 
cess using the vibrational state basis of the three elec- 
tronic potentials, the treatment would involve compli- 
cated transfer processes between a large number of states. 
Our presentation shows that the process can be under- 
stood very well in terms of the light-induced potentials, 
and their time-dependent vibrational eigenstates. Fur- 
thermore, this description also allows one to identify 
how the system parameters should be set in order to 
achieve specific outcomes. The main ingredient is estab- 
lishing the right balance between the adiabatic following 
of the time-dependent vibrational eigenstates, and nona- 
diabatic transfer between them at avoided crossings. We 
have discussed the situation in a molecular multistate 
framework, but in the light-induced potential description 
the relevant process really takes place within a general 
two-well potential structure. Thus our observations hold 
also for any two-well structure, where the barrier height 
and the well depths can be controlled time dependently 
in similar manner. 

This work was supported by the Academy of Finland, 
Project no. 43336. 
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